Potential of an ionic impurity in a large 4 He cluster 
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This paper presents an analysis of the motion of an impurity ion in a nanometer scale 4 He cluster. 
Due to induction forces, ions are strongly localized near the center of the cluster, with a root 
mean squared thermal displacements of only a few A. The trapping potential is found to be nearly 
harmonic, with a frequency of 2.3(1.0) GHz for a positive (negative) ion in a He cluster of radius 
5nm. The anharmonicity is small and positive (energy increases slightly faster than linear with 
quantum number). It is suggested that by using frequency sweep microwave radiation, it should be 
possible to drive the ion center of mass motion up to high quantum numbers, allowing the study of 
the critical velocity as a function of cluster size. 



I. INTRODUCTION 

The last few years have seen dramatic advances in the 
spectroscopy of atoms and molecules attached to large He 
clusters Q| . These clusters provide a unique environment 
for a spectroscopy which combines many of the attractive 
features of both high resolution gas phase spectroscopy 
and traditional matrix spectroscopy ||]. These include 
the ability to obtain rotationally resolved spectra of even 
very large molecules ||, and the ability to form and sta- 
bilize extremely fragile species [|J, including high spin 
states of molecules. 

Despite rapid progress, many fundamental questions 
remain about spectroscopy in this environment. One 
important topic that has received little attention is the 
dynamics of the center of mass motion of an impurity 
in a 4 He cluster 0. A recent experiment reported by 
Grebenev, Toennies, and Vilesov |(| has convincingly 
demonstrated that these clusters are superfluid. Thus 
we expect that an impurity, like in bulk superfluid He, 
should be able to move with little or no friction as long 
as its velocity stays below a 'critical velocity' which is 
found to be on the order of 30 m/s for motion of positive 
ions in bulk superfluid 4 He |?]||. Doped He clusters pro- 
vide an attractive system to study the size dependence 
of superfluid hydrodynamics. 

There have been several theoretical studies reported 
that consider the motion of an electron bound to the 
surface of He or other dielectric clusters ]9|-|lH . Most of 
these have only considered the electrostatic potential for 
the ion outside of the cluster, and solved for the quantum 
levels of the electron, including determining the minimum 
size for a cluster to bound an electron. The only excep- 
tions in the open literature known to the authors is the 
work of Antoniewicz et al. Jl(],[llJ that presents the elec- 
trostatic potential inside the cluster, but the potential 
given in that work is in error by a factor of two, and a 
paper by Northby, Kim and Jiang ]l5| that presents an 
approximate potential, similar to that given in section 2 
below, but without any derivation. The exact electro- 



static expressions have been given in the Ph.D. thesis of 
Ballester @ and Kim @. 

This paper will present a realistic potential for the mo- 
tion of an ion in a He cluster, based upon the electrostatic 
potential produced by the dielectric response of the He 
to the ion. Similar response is partly responsible for the 
trapping potential of an ion beneath a He-vapor inter- 
face that was exploited by Poitrenaud and Williams 
to determine the effective mass of positive and negative 
charge carriers in bulk He. This paper will then propose 
experiments that can be viewed as the natural extension 
of this earlier work on bulk He. In addition to allowing 
the determination of the size dependence of the effective 
mass of ions, it should also allow for the study of the 
effective size dependence of critical velocity of ion in the 
cluster. 

In a following paper, the motion of a neutral impu- 
rity atom or molecule will be considered, and an effective 
Hamiltonian for its motion, including both long range po- 
tential and hydrodynamic contributions, will be derived. 
The potential localizes the neutral impurity near the cen- 
ter of the cluster, though much more weakly than for the 
case of an ion which has much stronger long range in- 
teractions with the He, falling as r -4 versus r~ 6 . In the 
case of molecules, there is a coupling of the rotation and 
the center of mass motion which leads to a broadening 
mechanism in the rotational or ro-vibrational spectrum. 

II. POTENTIAL FOR AN ION IN A HE 
DROPLET 

A large fraction of the studies of impurities in He clus- 
ters have exploited mass spectrometry ||r9|-|2l||. The mi- 
gration of the charge in a cluster (which is most likely 
initially localized on a He atom) will be influenced by 
the effective potential of such a charge p2] |. By electro- 
statics, the charge will be most stable at the center of 
the cluster, furthest from the polarization charge that 
will develop on the cluster surface, due to the dielectric 
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constant of He. The purpose of this section is to derive 
an expression for this potential and explore some of its 
predicted consequences. This potential does not include 
the energy of solvation of the He around the impurity. 

Assume that we have an ion of unit charge, e, at a ra- 
dius a from the center of droplet, and pick the coordinate 
system such that the z axis is along the displacement of 
the ion from the center. In order to calculate the en- 
ergy, we can sum up the ion induced dipole interaction 
of the ion with the 'missing' He that would be 'outside' 
the droplet. This makes the 'zero' of energy an ion in an 
infinite bulk of liquid He, and also avoids the difficulties 
with the nearby He atoms which are strongly bound to 
the ion. 

Let r(9) be the distance from the ion to the droplet 
surface at a polar angle 9 measured from the ion. Basic 
trigonometry gives: 

R 2 = (a + r cos(6>)) 2 + r 2 sin(6>) 2 = a 2 + r 2 + 2ra cos(6) 

(1) 



from which we can derive: 



r(6) = \J R 2 — a? sin(6>) 2 - acos{6) 
The field from a charge is given by: 

E(r) 
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This leads to an energy difference from the bulk: 
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where p is the number density of He (0.0218 A |23|) 
and a is the polarizability of He. Evaluating the integrals 
(which were done using the Mathcad program [Q), gives 
the following result: 
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As the ion approaches the surface of the cluster (i. e. y — > 
1), the first term in the expression for Fo(y) dominates. It 
is easily shown that for e r ~ 1 and as y — > 1, the potential 
becomes the same as the 'image charge' potential for an 
ion approaching a planar He-vacuum interface from the 
Helium side [@. For liquid He at 3 K, e r = 1.05646 |§6). 
If we take R = 5 nm (which corresponds to ~ 11, 500 He 
atoms), the prefactor in the equation for AE/hc = 62 
cm -1 . If we compare the energy of an ion at the center 
of the cluster, y = 0, with one near the edge y = 0.90, we 
get a wavenumber difference of 151 cm -1 , which is very 
large compared to ki,T c /hc = 0.26 cm -1 for a 4 He clus- 
ter at T c = 0.38 K. Spectroscopic studies of a number 
of impurities have demonstrated that He clusters main- 
tain themselves, by evaporation, at a temperature close 
to this value |^7],[| , as had previously been predicted [^3| . 
This potential means that an ion will be strongly pushed 
towards the center of the cluster. Once in the center of 
the cluster, the thermal motion should produce a Gaus- 
sian distribution with RMS displacement of only 0.5 nm. 
This can be compared to the size of the 'snowball' of 
frozen He around a positive ion, which is known to have 
a radius of w 0.6 nm [p8[ . Electrons in He form a 'bub- 
ble' with a radius of « 1.7 nm p8fl . The mechanism and 
thus time scale for equilibration of an impurity center of 
mass motion with the internal motion of the He cluster 
is presently unknown, but most likely involves exchange 
of energy and angular momentum with quantized sur- 
face capillary waves, known as ripplons, as these are the 
only He cluster modes thermally excited at this temper- 
ature | p3fl . 

This treatment of the effective potential for an ionic 
impurity leaves out the interaction between the induced 
moments created by the charge. As long as one works 
in linear response theory, and one can treat the He as a 
continuum, then these effects can be included by a clas- 
sical electrostatic calculation, which are presented in an 
appendix. Because the relative dielectric constant, e r , is 
so close to one for liquid He, this more exact treatment is 
in excellent agreement with the more approximate treat- 
ment given above. This supports not only the neglect of 
three body effects for the present case of an ionic impu- 
rity, but also for the neutral impurities to be considered 
in a latter paper. 



F {y) ^l + \y 2 + \y A + 
3 5 



(7) 



In this paper, we will consistently use y for the reduced or 
fractional radius of the impurity ion or molecule. Using 
the relationships between the polarizablity, the electric 
susceptibility x, and the relative dielectric constant e r : 
pa = e x for \ < 1, \ = e r - 1, and e = e r e we 
can write this as: 
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A. Proposed Experiments 

Let us consider an ion in a He cluster of R = 5. Taylor 
expansion of the potential (Eq. |l^) around the center of 
the cluster gives a harmonic force constant of 
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If we assume an effective mass of 45 times the mass 
of 4 He, which was found in bulk He for positive ion 
mobility Jlq] , we get an effective vibrational frequency 
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v = 2.3 GHz. The dimensionless length corresponding to 
this vibration is 1.6 A. The vibrational transition mo- 
ment for the n= — > 1 transition is /iqi = 5.2Debye. 
The zero point level corresponds to a root mean squared 
(RMS) velocity of 3.8 ms -1 . It will require 90 quanta of 
vibration to reach a velocity of 30m s _1 , which is near the 
critical value for the onset of dissipation in bulk super- 
fluid He j?],U ■ Vibration around the center has a positive 
anharmonicity with x = 6.4 • 10~ 4 . This implies that 
the 89 — ► 90 vibrational transition will be blue shifted 
by fss 11%. Even for this high level of excitation, the 
classical turning point for the vibration is s=s 1.5 nm, and 
thus the ion will remain localized close to the center of 
the cluster, and thus away from the surface, where the 
approximations made in this paper are expected to break 
down. 

This suggests the following experiment to measure the 
cluster size dependence of the critical velocity. For ion 
motion with a peak classical value below the critical ve- 
locity, one expects little damping of the vibrational mo- 
tion. Thus, excitation of the vibrational transition should 
have little or no observable effects. However, once the vi- 
brational velocity of the ion exceeds the critical velocity, 
the motion will become strongly damped and the cluster 
will continue to absorb energy from a microwave held. 
This will lead to evaporation of He atoms which can be 
detected by a mass selected and mass analyzed beam. By 
chirping the microwave frequency, one can exploit the an- 
harmonicity to produce an almost pure number state of 
the motion, at least before damping becomes important, 
with the level of excitation determined by end point of 
the chirp. It is easily verified that due to the large tran- 
sition dipole moments, only modest microwave power is 
required to drive the ion up to high levels of excitation. 
The motion of positive ions in bulk He is associated with 
vortex rings, which have the curious property that their 
velocity is inversely proportional to their energy |t| . The 
effect of confining the ion in a nanoscale He cluster could 
dramatically alter the dynamics of ion motion. 

Even below the critical velocity, one may expect some 
coupling between the motion of the ionic impurity and 
the internal degrees of freedom of the He cluster, par- 
ticularly the surface ripplons, which are low frequency 
motions. Using the methods described in the paper that 
follows this one, on the potential of a neutral impurity 
in a large He cluster, it is possible to estimate the size of 
the coupling of the ionic motion to the ripplons. These 
couplings will be expected to cause perturbations of the 
frequency of center of mass motion of the ion as a function 
of excitation level and/or cluster size, as the ion motion 
and ripplons pass through resonance conditions. It is eas- 
ily seen that the lowest order in the interaction energy 
between the ion displacement, a, and a displacement of 
a ripplon of angular momentum L, S(L), is proportional 
to: 

H^^^) e -^(a/R^S(L) (10) 



This and higher order coupling terms will allow energy to 
flow from the ion motion to the ripplons. Energy in the 
ripplons, in turn, can lead to evaporation and observable 
reduction in the size of a mass selected cluster. Such 
experiments could provide the first measurements of the 
excitation levels of nanometer scale He clusters. 

Another potential experiment is to examine the reso- 
nances of electrons in He clusters. Such solvated elec- 
trons are metastable, and the electrons will eventually 
be expelled from the cluster [p9| . Kim, Yurgenson and 
Northby have observed by spontaneous and infrared 
induced electron detachments from large He clusters sev- 
eral milliseconds after formation. The stability of such 
charged clusters to electric fields demonstrates that the 
electrons are in bubble states of the He cluster, not at- 
tached as surface states pl[ . Excitation of the center 
of mass motion of these 'bubble states' should result in a 
dramatic increase in the rate of electron evaporation from 
the cluster. The potential for the negatively charged bub- 
ble should be the same as for a positive ion, as long as the 
center of the bubble is further from the surface than its 
radius, f» 1.7nm. The primary difference of the bubble 
compared to the positive ion is that the effective mass 
is w 243 times the mass of 4 He |3^| , which means that 
for a 5 nm cluster, the harmonic vibrational frequency 
is ss 1 GHz. Since the Harmonic vibrational frequency 
is proportional to the inverse square root of the cluster 
radius, it should be possible to selectively neutralize all 
clusters below a certain size by sweeping the microwave 
field from a certain frequency to higher values. By us- 
ing microwave double resonance experiments, it should 
be possible to determine the homogeneous width of the 
resonance, and thus any drag that the electron bubble 
inside the cluster may experience. 

III. SUMMARY 

This paper has developed the potential that governs 
the motion of an ion inside a nm scale He cluster. This 
motion is determined by long range electrostatic interac- 
tions, and thus can be calculated without having to deal 
with the much more difficult problem of the energetics of 
solvation of the ion in liquid He. It is found that despite 
the small dielectric constant of He, the potential is very 
effective at driving the impurity ion into the center of 
the cluster. In thermal equilibrium, RMS displacements 
of only a few A are expected. 

After ionization by electron impact, the He ion under- 
goes a 'random walk' due to charge exchange with other 
He atoms until it becomes localized as a HeJ core which 
is then solvated (3^] . The electrostatic potential derived 
here will strongly bias this 'random walk', given the low 
temperatures of the clusters j22|. Such a bias could be 
expected to effect the dynamics following ionization, in- 
cluding the probability of charge exchange with an im- 
purity, which will also likely be localized near the center 
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of the cluster. 

The trapping potential of an ion is found to be highly 
harmonic, with a frequency in the low microwave region, 
and having very large transition dipole moments. Exper- 
iments are suggested that would exploit these resonances, 
and the ability to drive the ion up to high quantum states 
by an adiabatically sweep microwave pulse. 
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We can determine A n and B n by satisfying the bound- 
ary conditions V (R,9) = Vi(R,0) and e E ^ r — eEi^ r . 
For these equations to hold for all 6, the respective co- 
efficients of each P n must be equal. This leads to the 
following expressions: 
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We can now calculate the electrostatic energy by using 
the expression given in Eq. 4.83 of Jackson pa): 



APPENDIX: EXACT ELECTROSTATIC 
TREATMENT FOR ION IN HE CLUSTER 

In this appendix, the classical electrostatic energy that 
determines the motion of an ion inside a He cluster will be 
derived, modeling the He cluster as a dielectric sphere. 
Given the errors in the previous open literature, it is 
important that the expressions given below be justified. 
Further, comparison of the exact and approximate energy 
expression given by Eq. || is important since it allows for 
an estimate to be made of the approximations used in 
the companion paper for the energy of neutral impurities 
in He, where a treatment including many body effects is 
not practical. 

Outside of the cluster, we must satisfy Laplaces equa- 
tion and we must have only terms which go to zero at 
infinity. Thus: 



V (r,6) = J2 B nr- {n+1) Pn(cost 



(11) 
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Inside the cluster, we have a field due to a point charge 
(screened by the dielectric) plus fields due to surface 
charges. The field of the later must also satisfy Laplaces 
equation, but with positive powers: 



Vi{r,0) 



e 1 

47re \r — a\ 



J2A n r n P n (cos6) (12) 



n=0 



to determine the angular dependence of the first term, 
we use the expansion: 
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where r< = min(r,a) and r> = max(r,a). We will be 
primarily interested in the potential for r > a 
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where p(x) = ed(x — a) is the free charge density. Since 
our reference energy is that of a point charge in an infinite 
dielectric, we must subtract off the monopole term in the 
potential. This gives: 
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Putting in the above equation for A n we get: 
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This potential is exactl y o ne half that previously reported 
by Antoniewicz et al. [Tcjjll]]. These authors reference a 
classic text [ [34) for their expression, but unfortunately, 
the present authors could not locate a copy of this book. 

In the limit that e r —> 1, Eq. [l9| agrees exactly with the 
Taylor expansion of the expression given above (Eq. ||) 
for the energy obtained by 'adding up' the ion-induced 
dipole energy contributions. Since this earlier expression 
correctly gives the correct planer 'image charge' poten- 
tial as y — * 1, it cannot be in error by a factor of two, 
and thus the error must lie with the expression given 
by Antoniewicz, Ballester et al. The Ph.D. thesis of 
Ballester [jj"6f gives the correct power series expression 
without the factor of two error. The correct expression 
was also given in the Ph.D. thesis of Kim Jl7[] , who cor- 
rected some other minor errors in Ballester's expressions. 

For the case of He clusters, if the final sum is dropped 
in Eq. Q the resulting approximation is in error by at 
most few percent over the range < y < 0.9. In this same 
range of y, Eq. || agrees with Eq. 19 even better, with 
error of less than 1%. Thus the fractional error of the 
approximate treatment that neglects many body effects 
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is considerably smaller than might have been predicted 
a priori, w (e r — 1). 

The power series expansion for W (Eq. |l9| ) is com- 
putationally useful for small y values (say < 0.5), but 
convergence slows dramatically as y — > 1, since the real 
solution diverges in that limit. By reference to the form 
of the correct solution in the limit of e r » 1 (Eq. ||), it is 
possible to subtract out the power series of the divergent 
parts of the solution. The resulting expression is: 
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In this form, the fact that W goes exactly into Eq. |B| in 
the limit e r — > 1 is transparent. Further, the last sum is 
convergent, even at y = 1, though the convergence is slow 
there. However, in the region of slow convergence, the 
sum makes a negligible contribution to the total energy. 
Further, if the sum is truncated at the n — N — 1 term, 
the remainder of the sum can be approximated by an 
integral: 
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(22) 



where T{a,x) is the incomplete Gamma function. 
Ballester |l6| suggests using the method known as 
Aitkin's S 2 process J3f| to speed up convergence of the 
series. Kim Jl7| derived a similar closed form expression 
for the interaction energy by a different method. 
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